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SUMMARY 

A simplified model is analyzed in order to give some insight into the effects of buoy- 
ancy and shear flow on turbulence. Two- point correlation equations, which contain mean 
velocity and temperature gradients, as well as body force terms, can be constructed 
from the Navier-Stokes, energy, and continuity equations. While previous papers by the 
author considered the effects of shear and buoyancy separately, the present paper con- 
siders their combined effects. In that case the ratio of buoyancy to shear effects, as 
given by the Richardson number, is a consideration. The velocity and temperature gra- 
dients, as well as the body force, are considered to be vertical and uniform. 

An initial value problem in which the turbulence is initially isotropic is solved to il- 
lustrate the effect of buoyancy and shear parameters on the turbulence. The initially iso- 
tropic turbulence quickly becomes anisotropic under the influence of the buoyancy and 
shear. Although the turbulence is assumed to be weak enough for terms in the equations 
containing triple correlations to be negligible, a term proportional to the velocity gradi- 
ent transfers energy between wave numbers. Various components of the mean turbulent 
fluctuations and their spectra, as well as the ratio of eddy conductivity to eddy viscosity, 
are calculated as functions of buoyance and shear parameters. For highly stabilizing 
conditions, the interaction of the shear and buoyancy forces with the turbulence can pro- 
duce negative eddy conductivities and viscosities. The analytical results are compared 
with available experimental data. 


INTRODUCTION 

The effects of buoyancy and of shear on weak homogeneous turbulence are considered 
separately in reference 1 and 2. In real situations, for instance in the atomsphere, the 


two effects often occur simultaneously. The speculative theories given in references 3 
and 4 consider that case. 

In the present paper, the methods used in references 1 and 2 are extended to analyze 
the combined effects of buoyancy and shear on homogeneous turbulence. Correlation equa 
tions for velocities and temperatures at two points in the turbulent field are constructed 
from the Navier- Stokes, energy, and continuity equations. In order to obtain a determi- 
nate set of equations, the turbulence is assumed to be weak enough to neglect the terms 
in the equations containing triple correlations. Although this assumption may limit the 
upper Reynolds number for which the analysis is valid, the analysis is of interest in itself 
in that it gives an asymptotically exact solution for turbulence in the limit of low Reynolds 
numbers. Moreover, when mean velocity and temperature gradient terms are present in 
the equations, the turbulence may not have to be as weak as it would otherwise for the 
triple correlation terms to be negligible compared with the other terms. As will be seen 
later, comparison of the analysis with experiment indicates that the theoretical results 
bear a correspondence with real turbulence at moderate Reynolds numbers. 

The correlation and spectral equations required in the analysis will be considered in 
the next section. 


BASIC EQUATIONS 


The Navier- Stokes equation with buoyancy effects included can be written as (ref. 5) 
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where the subscripts (except e) can take on the values 1, 2, or 3, and a repeated sub- 
script in a term signifies a summation. The quantity u ^ is an instantaneous velocity 
component, T is the instantaneous temperature, x^ is a space coordinate, t is the time, 
p is the density, v is the kinematic viscosity, p is the instantaneous pressure, g^ is a 

component of the body force, and 0 = -(l/p)(0p/9T)_ is the thermal expansion coefficient 

• P 

of the fluid. The quantities T g and p g are the equilibrium temperature and pressure, 
respectively. (All symbols are given in the appendix. ) In obtaining the last term in equa- 
tion (1) (buoyancy term), the density is assumed to depend effectively only on tempera- 
ture and is not far removed from its equilibrium value (value it would have for no heat 
transfer). Note that the equilibrium temperature is uniform whereas the equilibrium 
pressure is not. The instantaneous velocities and temperatures in equation (1) can be 
divided into mean and fluctuating components as follows: 
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Then, if equation (1) is averaged, and the averaged equation is subtracted from the unav- 
eraged one, 
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Taking the divergence of equation (4) (differentiating with respect to x.) and using continu- 
ity give 
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The instantaneous energy equation can be written as 
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Substituting equations (2) and (3) in equation (6), averaging, and subtracting the averaged 
equation from the unaveraged one give 
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Equations (4), (5), and (7) apply at a point P in the turbulent fluid. Similar equations at 
another point P’ can be obtained simply by priming the variables and changing the sub- 
script i to, say, j. Equations involving correlations between fluctuating quantities at 
P and P’ can then be constructed by methods similar to those used in references 1, 2, 
and 6. In fact, the equations considered in those references are special cases of those 
obtained here. The resulting equations for homogeneous turbulence with uniform velocity 
and temperature gradients are 
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Equations (8) to (15) form a determinate set if the turbulence is weak enough to ne- 
glect terms containing triple correlations in comparison with the other terms. They can 
be converted to spectral form by introducing three-dimensional Fourier transforms de- 
fined as follows: 
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where IT is the wave number vector and d/T = dic^ d d k^. The magnitude of k has the 
dimension 1/length and can be considered to be the reciprocal of a wavelength or eddy 
size. Then, from equation (16) (see eq. (10), ref. 7), 


5 


(24) 


5 1K • r dr 


The remainder of the paper will present the case where the velocity and the temper- 
ature gradients are in the Xg- direction (vertical) and the body force (gravity) is in the 
-Xg-direction. Then, let 

g - -g 3 (2 


and neglect terms containing triple correlations; equations (8) to (15) become, in spectral 
form, 
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where 6 - is the Kronecker delta. Substitution of equations (32) to (35) into equa- 
tions (28) to (31) shows that cp^ = cp^ and y. = y! for all times if they are equal at an 
initial time. It will be assumed herein that the turbulence is initially isotropic and that 
the temperature fluctuations are initially zero, so that the above relations will hold. 
Thus, the set of equations (28) to (35) becomes 
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Equations (36), (37), and (38) give contributions of various processes to the rates of 
change of spectral components of mm, ruT, and r2, respectively. The second term in 
each equation is a transfer term which transfers activity into or out of a spectral compo- 
nent by the stretching or compressing of turbulent vortex filaments by the mean velocity 
gradient, as discussed in references 2, 6, and 7. The terms with k in the denomina- 
tor are spectral components of pressure-velocity or pressure-temperature correlations 
and transfer activity between directional components (ref. 2). The terms proportional 
to /3g and 6 - 3 (or 6 j 3 ) are buoyancy terms which augment or diminish the activity in a 
spectral component by buoyant action. The last terms in the equations are dissipation 
terms, which dissipate activity by viscous effects (eq. (36)) or by conduction effects 
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(eq. (38)). The dissipation term in equation (37) contains both viscous and conduction ef- 
fects, inasmuch as it dissipates spectral components of velocity-temperature correla- 
tions. The remaining terms in the equations produce activity by velocity or temperature 
gradient effects. Although a buoyancy term does not appear in equation (38), buoyancy 
affects 6 (or t2) indirectly through the temperature gradient and (or rug). 

For solving equations (36) to (38), the turbulence is assumed to be initially isotropic 
at t = tg. That condition is satisfied by the relation 
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where Jq is a constant that depends on initial conditions (ref. 2, eq. (43)). For the ini- 
tial conditions of 6 and y^ (at t = tg), it is assumed that 

5 0 =(yi) =0 (40) 
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That is, the turbulence- producing grid is assumed to be unheated, so that the temperature 
fluctuations are produced by the interactions of the mean temperature gradient with the 
turbulence. 


SOLUTION OF SPECTRAL EQUATIONS 


In order to reduce the set of partial differential equations, equations (36) to (38), to 
ordinary differential equations, the running variables £ and t\ are considered, of which 
Kg and t are particular values such that £ = Kg when rj = t. If | and r\ are intro- 
duced into the set of equations in place of Kg and t, the resulting equations will, of 
course, automatically satisfy the original set. In addition, 
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where the subscripts outside the parentheses signify quantities which are held constant. 
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Then equations (36) to (38) become 
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Note that the first three of these equations are independent of the remaining ones. 

The constant in equation (41) subject to the initial conditions, equations (39) and (40), 
may be determined by letting £ = when V ~1q> or £q = £ + a Kj(rj - rj 0 ). This equa- 
tion applies for any value of £ , and thus 
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Equation (50) gives the value of t at which to start the integration for given values of k 3 , 
K a, and t-t Q . The initial conditions, equations (39) and (40), may be satisfied by letting 
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The integration of equation (43) to (49) then goes from to £ = Kg. Final values of 
( p i y., and 6, for which £ = k 3 and r\ = t are of interest. The quantity t, can be con- 
sidered as a dummy variable of integration. 

The following dimensionless quantities are introduced in order to convert equa- 
tions (43) to (49) to dimensionless form: 
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In addition, spherical coordinates are introduced into the equations by using the trans- 
formations 
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Equations (43) to (49) then become 
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For integrating these equations, £* starts at 


= k* cos 0 + a*K* cos cp sin 0 


where 
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The integrations were carried out numerically on a high-speed computer for various fixed 
values of k*, 9, <p, a, and Ri. Directionally integrated spectrum functions can be ob- 
tained from (refs. 6 and 7) 



In this equation, 


is the vorticity spectrum tensor given by reference 8. 
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The primed quantities <p'^ and S2!. give, respectively, components of cp^ and fi^ in 
a coordinate system rotated 45° about the x 2 -axis from the x^-axis toward the Xg-axis. 
Since cp .. and fi^ are second order tensors, components in the rotated system are 
(ref. 9) 
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The spectrum functions given by equation (67) can be integrated over all wave num- 
bers to give 



(71) 


where the primes again refer to components in a coordinate system rotated 45°. Com- 
puted spectra and correlations will be considered in the next section. 


RESULTS AND DISCUSSION 


All the calculated results given in this section are for a gas with a Prandtl number 
of 0.7. Dimensionless energy spectra ^spectra of UjU|j and spectra of t 2 are plotted 
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in figures 1 and 2. The spectra are plotted for several values of a* = (t - tQ)dUj/dx 3 
and of g* = /3g(t - tp) 2 dT/dx^. The parameter g* rather than the Richardson number 
Ri = ;3g dT/dXg /(dUj/dXg)^ is used here since the use of g* and a* enables us to con- 
sider buoyancy and shear effects separately. (The Richardson number contains both 
buoyancy and shear effects. ) The quantity g* is related to Ri and a* by the equation 

g* = a* 2 Ri (72) 

When plotted by using the similarity variables shown in figures 1 and 2, the dimen- 
sionless spectra for no buoyancy and shear effects g* = a* =0 do not change with time, 
and thus comparison of the various curves indicates how buoyancy and shear effects will 
alter the spectra. Thus, if a dimensionless- spectrum curve lies above the curve for 
g* = a* = 0, the turbulent activity for that case is greater than it would be for no buoy- 
ancy and shear effects. The turbulence itself decays with time, as in references 1 and 2. 

Positive values of the buoyancy parameter g* correspond to stabilizing conditions, 
and negative values correspond to destabilizing conditions. Figures 1 and 2 show that the 
trends with g* for a case with shear (a* = 2) are similar to those from reference 1 for 
no shear (a* = 0). That is, in the destabilizing case, buoyancy forces tend to feed energy 



Fi gure 1. - Dimensionless spectra of turbulent energy 
UjUj*. Prandtl number, 0.7. 



Figure 2. - Spectra of dimensionless temperature vari- 
ance -ft . Prandtl number, 0.7. 
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or activity into the turbulent field, whereas in the stabilizing case they tend to extract it. 
Comparison of the curves with shear (solid curves) with those without shear (dashed 
curves) for values of g* of -4, 0, and 4 indicates that for all three cases the effect of 
the shear is to feed energy or activity into the turbulent field. Thus for the destabilizing 
case, the buoyancy and shear have similar effects; but for the stabilizing case, they have 
opposite effects. Comparison of the curve in figure 1 for g* = 4 and a* = 2 with that 
for a* = g* = 0 indicates that for the former curve the energy added by the shear effects 
approximately balances that extracted by buoyancy but that the wave number distributions 
for the two processes are slightly different. 

As a* increases, the spectra become asymmetric, the slopes on the high-wave- 
number sides of the curves becoming more gradual. The dot- dashed curves for an a* 
of 4 and a g* of 1.6 are plotted to show this effect. As in references 2, 6, and 7, the 
effect is due to the transfer of energy or activity into the high- wave- number regions by 
the transfer term associated with the mean velocity gradient (see the discussion in the 
paragraph following eq. (38)). 

The buoyancy forces might be expected to act more strongly on the turbulent velocity 


components lying in the direction of those forces than on the other components. This ex- 
pectation is confirmed in the plot of u j|/ u ^ for a* = 0 in figure 3. The ratio u |/ u j is 
greater than 1 in the destabilizing case and less than 1 in the stabilizing case. However, 
although u|/ u| becomes small, it does not approach zero for highly stable conditions. 


Apparently begins to decrease as rapidly as, or more rapidly than u|, as g* be- 
comes large. __ 

The voriticity component ratio u)f|/ is also plotted in figure 3. The trends for 

0,2/cof are opposite to those for u|/uj. For 



Figure 3. - Velocity and vorticity-variance component 
ratios plotted against buoyancy parameter for case of 
no shear. Prandt I number, 0.7. 


the stabilizing case, the voriticity tends to be 
alined in the direction of the buoyancy forces. 

The turbulent velocities associated with that vor- 
ticity will then tend to be normal to the buoyancy 
forces, in agreement with the curve for u^/u^. 
For the destabilizing case, the vorticity tends 
to lie in directions normal to the buoyancy forces. 
That will tend to increase u|/ u| as shown in 
figure 3, although the ratio will not approach in- 
finity, because even if all the vorticity lies in 
directions normal to the buoyancy forces, u | or 
u| will not be zero. 

In the preceding discussion, the effect of 
buoyancy forces were considered on the turbu- 
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Figure 4. - Velocity- and vorticity-variance com- 
ponent ratios in coordinate system rotated 45° 
plotted against shear parameter for case of no 
buoyancy effects. 


lence components without mean shear. Next, 
consider the case of shear with no buoyancy 
effects. In that case, the turbulent vorticity 
(or vortex filaments) would be expected to 
tend to aline in the direction of maximum 
strain, which is at 45° to the mean velocity. 
Figure 4 shows turbulent vorticity and velocity 
components ratios in a coordinate system ro- 
tated 45° counterclockwise about the X 2 ~axis. 

If the vorticity were all alined in the direction 
of maximum strain, cof^/co^’ would be zero. 
The curve shows that there is a strong tend- 
ency for that alinement to occur at moderate 
values of a*, but the degree of alinement 
does not continue to improve as a* becomes 
large. The tendency for the vortex filaments 
to aline in the direction of maximum strain is 
reflected in the trend for the turbulent velocity 
components to become maximum in a direc- 
tion normal to the maximum strain, as shown 
in the curve for u|/ u|\ The degree of 
alinement, however, does not continue to im- 
prove as a* becomes large. 

Combined effects of buoyancy and shear 
on uj^/ u| are shown in figure 5. The curves 
show that for no buoyancy effects (g* = 0) the 
turbulence component u|, which is in the di- 
rection of the mean velocity gradient, is re- 
duced in comparison with u^ by the shear. 
This trend also occurs for negative (destabil- 
izing) values of g* and for small positive 

(stabilizing) values of g*. For more strongly stabilizing conditions, the trends become 
more complex, and the curves cross over one another. 

The effects of buoyancy and shear are considered separately in figure 5, which uti- 
lizes the parameters g* and a*. Since the Richardson number contains both buoyancy 
and shear effects, one might suppose that its use would reduce or eliminate the need for 
another parameter. Figure 6 shows that is not the case, since u^/ u^ is a strong func- 
tion of both Richardson number and a*. 

The ratio of two turbulence components which are normal to the body forces and 
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Figure Plot showing velocity-variance component 

ratio u^/u| as a function of buoyancy and shear 
parameters. Prandtl number, 0.7. 
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Figure Plot showing velocity-variance component 
ratio uyu| as a function ot Richardson number 
and shear parameter. Prandtl number, 0.7. 



Figure 7^_- Plot showing velocity-variance component 

ratio u^/uj as a function of buoyancy and shear 
parameters. 


mean velocity gradient plotted as a function of g* and a* are shown in figure 7. For 
no shear (a* = 0), u|/ is 1 since the turbulence is axially symmetric. For a* * 0, 
the shear tends to destroy the axial symmetry and to reduce xs^j u^ below 1. 

Consider next the turbulent heat transfer and the turbulent shear stress. 

//“ 2 \ 1 / 2 /“ 2\ 1 / 2 

Temperature-velocity correlation coefficients -TUg/lr I l Ug I are plotted in fig- 
ure 8. The correlation TUg is proportional to the turbulent heat transfer in the direction 
of the temperature gradient. The unusual feature of these results is that TUg changes 
sign as g* becomes large. That is, for very stable conditions, the turbulence begins to 
pump heat against the temperature gradient. This phenomenon was observed in the re- 
sults of reference 1 for an a* of zero. As a* increases, the value of g* at which 
TUg changes sign increases. 
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Figure 8 . - Plot showing temperature-velocity correlation 

coefficient -iuj/t^u?^ as a function of buoyancy 
and shear parameters. Trandtl number, 0.7. 
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Figure 9._^Piqt_$howing velocity-velocity correlation coefficient 
-u^/ufup as a function of buoyancy and shear para- 
meters. 


Velocity-velocity correlation coefficients for shear are plotted in figure 9. At small 
values of g*, the trends with a* are opposite to those for figure 8; the values of 
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are 


close to 1 for small a* and g*. As was the case for TUg, u^Ug changes sign as g* be- 
comes large. As conditions become strongly stabilizing, the turbulence begins to pump 
the fluid in such a way as to tend to increase the velocity gradient. Thus, there occurs, 
for sufficiently large values of g*, a negative eddy viscosity as well as a negative eddy 
conductivity. Although a negative eddy conductivity can occur with only buoyancy effects 
present, the occurrence of a negative eddy viscosity requires combined buoyancy and 
shear effects. 

A possible explanation for the theoretically observed negative eddy viscosity and con- 
ductivity can be given in terms of a modified mixing length theory as illustrated in 
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Distance, 


figure 10. Normal turbulent heat transfer 
or shear stress is shown in the left portion 
of the figure. An eddy originating at the 
mean velocity and temperature of the fluid 
at a point may move either upward or down- 
ward. By conduction and viscous effects, it 
will tend to acquire the local mean temper- 
ature and velocity of the fluid as it moves, 
and thus its path will curve toward the mean 
velocity or temperature line. When the eddy 
mixes with the fluid, it will tend to decrease 
the mean temperature and velocity gradients. 
The effective eddy conductivity and viscosity 
will be positive since they act in the same 
direction as the molecular conductivity and 
viscosity. 

By contrast, for the abnormal case 
where the buoyancy forces are strongly stabilizing, the original direction of motion of an 
eddy may be reversed. This reversal might happen because the buoyancy force, in the 
stabilizing case, acts in the direction opposite to that in which the eddy starts to move. 
Possible paths for the eddy on the distance-temperature or distance-velocity plane under 
these conditions are sketched on the right side of figure 10. As shown, the eddy path can 
cross the mean temperature or velocity line. As the eddy mixes with the fluid, it will 
then tend to increase the mean temperature and velocity gradients, and thus the effective 
eddy conductivity and viscosity will be negative. The actual mechanism may be more 
complicated than that considered here. The preceding explanation is given only to show 
that negative eddy conductivities and viscosities are physically reasonable. The turbulent 
heat transfer and shear stress do not necessarily change sign at the same value of g*, 
since the eddy paths on the distance-temperature plane and on the distance-velocity plane 
may be different because of differences between the conduction and viscous effects on the 
eddy as it moves. Comparison of figures 8 and 9 shows that the turbulent shear stress 
changes sign first as g* increases. 

The ratio of eddy conductivity to eddy viscosity plotted against g* is shown in fig- 
ure 11. The eddy conductivity and eddy viscosity are defined by the relations 



Figure 10. - Sketch illustrating possible mechanism for pro- 
ducing negative eddy conductivity and viscosity. 
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and 


Figure 11. - Plot showing ratio of eddy conductivity 
to eddy viscosity as a function of buoyancy and 
shear parameters. Prandtl number, 0.7. 



Figure 12. - Plot showing correlation coefficient 
Yu^/ t 2 1,2 u 2 1 ' 2 as a function of buoyancy and 
shear parameters. Prandtl number, 0.7. 


dUj 

dx 3 

The ratio e ^/e is calculated from Tui^/^UjU^ /a*) . For small values of a*, e de- 
creases with increasing g* except for large g*. For an a* of 4, e^/e increases with 
increasing g*. The sharp increases in e h /e near the ends of the curves occur because 

the eddy viscosity approaches zero and changes sign near those points. 

//^\l/2/~2\l/2 

Values of the correlation coefficient TUj/l r ) UijJ are presented in figure 12. 
The correlation TUj is proportional to turbulent heat transfer in the x^-direction. The 
fact that there should be heat transfer in the Xj-direction is surprising since there is a 
temperature gradient only in the Xg-direction. It appears, however, that tuj can be 
nonzero because of the nonzero values of rug and UjUg. Since there is a correlation 
between t and Ug and between Ug and Uj, the fact that a correlation should occur be- 
tween r and Uj seems reasonable. It must be admitted, though, that heat transfer in a 
direction of zero temperature gradient runs contrary to normal intuition. It should be 
noted that the effect is not dependent on the presence of buoyancy forces (i. e. , g* can be 
zero). The turbulent heat transfer in the Xj-direction is not necessarily small compared 
with that in the Xg-direction. Figure 13, which shows plotted values of -TUg/ TUj, indi- 
cates that the turbulent heat transfer in the two directions can be of the same order of 
magnitude. 

Comparison of the present analytical results with available experimental data is of 
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Buoyancy parameter, g* - pg(t - t 0 ) 2 dT/dx 3 

Figure 13. - Plot showing temperature-velocity correlation 
component ratio -tUj/tUj as a function of buoyancy 
and shear parameters. 


interest in order to see if there is a corre- 
spondence. Experimental data for grid- 
generated turbulence in the presence of com- 
bined buoyancy and shear effects are pre- 
sented in references 10 and 11. A compari- 
son between analysis and the data obtained at 
station 5 in reference 11 is given in figure 14. 
The various turbulence ratios are plotted 
against Richardson number rather than 
against g*, since (t - t ) in g* is hard to 
estimate because of uncertainties in the initial 
and other conditions in the experiment. For 
the same reason, no attempt is made to cal- 
culate a* = (t - t ) dUj/dXg directly for the 
experiment. The data and analysis for 
u 3 / Ri = 0 indicate a value for a* of 

about 2, and that value is used throughout the 
comparison. A reasonable correspondence 



Richardson number, Ri 


Figure 14. - Comparison of analytical results with experimental data at 
station 5 in reference 11. Shear parameter, a* ■ 2. 
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exists between analysis and experiment for the various turbulence ratios and correlation 
coefficients except in the last plot. Conditions in the analysis and experiment differ. 

For instance, the analysis assumes that the temperature fluctuations arise from the in- 
teraction of the mean temperature gradient with the turbulence whereas in the experiment 
the heated grid probably produced temperature fluctuations. A disconcerting aspect of 
the experimental data is that it shows an increase in turbulence level with distance while 
the analysis shows a decrease. The increase shown in the experiment might be caused 
by lateral inhomogenieties in the turbulence, boundary layer effects, or by triple correla- 
tion effects. All these effects were neglected in the analysis. 


SUMMARY OF RESULTS 

The analytical results for combined effects of vertical buoyancy forces and vertical 
velocity gradients indicate that, as in the case of no shear, destabilizing buoyancy forces 
can feed energy or activity into a turbulent field whereas stabilizing buoyancy forces can 
extract it. The effect of the shear is to feed energy or activity into the turbulent field. 
Thus for the destabilizing case, the buoyancy and shear have similar effects; but for the 
stabilizing case, they work in opposite directions. 

Energy or activity transfer between wave numbers by the stretching of turbulent vor- 
tex filaments by the mean velocity gradient causes the spectra to become asymmetric; 
the slopes on the high- wave- number sides of the spectra become more gradual. 

For the destabilizing case, buoyancy forces tend to increase the vertical turbulence 
component in comparison to the horizontal component in the flow direction while the 
shear tends to decrease it. For weakly stabilizing conditions both the buoyancy and 
shear tend to decrease the ratio of vertical to horizontal turbulence components. For 
more strongly stabilizing conditions, the trends become less well defined. 

The shear tends to aline the turbulent vorticity in the direction of maximum mean 
strain, which is 45° from the flow direction. Destabilizing buoyancy forces tend to aline 
the vorticity in horizontal directions whereas stabilizing forces tend to aline it vertically. 

When buoyancy forces are strongly stabilizing, the eddy conductivity and viscosity 
can be negative. This result appears reasonable when considered from the standpoint of 
a modified mixing length theory. Turbulent heat transfer occurs in a horizontal as well 
as a vertical direction, even though the velocity and temperature gradients are both 
vertical. 

Comparison of analytical with experimental results indicates that the former bear a 
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reasonable correspondence with results for observed turbulence, although the conditions 
in the experiment differed somewhat from those assumed in the analysis. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 22, 1966, 

129-01-09-07-22. 
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APPENDIX - SYMBOLS 



a 

vertical velocity gradient, dUj/dXg 

y i 

defined by eq. (21) 

a* 

shear parameter, defined by 

y i 

defined by eq. (22) 


eq. (56) 

y* 

r i 

defined by eq. (54) 

b 

g 

vertical temperature gradient, 
dT/dx 3 

-gg, vertical body force/unit mass 

A 

6 

defined by eq. (67) 
defined by eq. (23) 


in -Xg- direction (gravitational 

6 « 

Kronecker delta 


force/unit mass) 

5* 

defined by eq. (55) 

g 3 

vertical body force/unit mass 

€ 

eddy viscosity 

g i 

body force component/unit mass 

e h 

eddy conductivity 

g* 

o 

buoyancy parameter, /3g(t - t Q ) 

c 

defined by eq. (19) 


dT/dx 3 

c* 

defined by eq. (20) 

J 0 

constant that depends on initial 
conditions 

V 

running or dummy variable that 
equals t when £ = 

Pr 

Prandtl number, v/a 

9,(p 

angular coordinates (see eq. (59)) 

P 

pressure 

K i 

wave number component 

F ’ r i 

vector between points P and P* 

K* 

defined by eq. (51) 

Ri 

Richardson number, defined by 

A ii 

defined by eq. (67) 

T 

t 

U i 

eq. (57) 

mean temperature 
time 

mean velocity component 

X i 

H 

V 

h 

defined by eq. (17) 
defined by eq. (18) 
kinematic viscosity 
running or dummy variable for 

u i 

fluctuating velocity component 


which Kg is a particular value 

iTuT* 

v 5/2 (t - t Q ) 5/2 

J Vi 

J o 

p 

density 

i J 

T 

temperature fluctuation 

x i 

space coordinate 


temperature-velocity correlation 

a 

thermal diffusivity 


- ,//2 aq 

P 

r j 

expansion coefficient 
defined by eq. (67) 

TU^ 1 

u i 

J 0 b 
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1 


in flow direction 


t 2* 

y 5/2 (t - 1 0 ) 1/2 

t 2 

1 

j 0 h 2 



defined by eq. 

(16) 

*8 

defined by eq. 

(53) 


defined by eq. 

(67) 

1 ] 

At - 1 0 ) 2 


*1] J 

J o 


n.. 

n 

given by eq. (68) 


co .ok turbulent vorticity variance 
Subscripts: 

0 at virtual origin of turbulence 

where turbulent energy would be 
infinite (It is assumed that turbu- 
lence is isotropic at Xq and 
that velocity and temperature 
gradients begin to act there) 


3 in vertical direction, which is di- 
rection of mean velocity gradient 
and buoyancy force 

Superscripts: 

* at point P', also refers to coordi- 

nate system rotated 45° 

* on dimensionless quantities 

— over averaged quantities 
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